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Abstract. We compare the magnetic field at the center of and the self-magnetic 
flux through a current-carrying circular loop, with those obtained for current-carrying 
polygons with the same perimeter. As the magnetic field diverges at the position 
of the wires, we compare the self-fluxes utilizing several regularization procedures. 
The calculation is best performed utilizing the vector potential, thus highlighting 
its usefulness in practical applications. Our analysis answers some of the intuition 
challenges students face when they encounter a related simple textbook example. These 
results can be applied directly to the determination of mutual inductances in a variety 
of situations. 
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1. Introduction 

A common exercise in introductory physics courses concerns the comparison between 
the magnetic fields due to two loops of equal length P, carrying the same current i, one 
shaped into a square and the other shaped into a circle. One is asked to compare the 
magnetic fields at the centers of the respective figures jl], finding that the field at the 
center of the square is larger than the field at the center of the circle. In our classes, 
this problem is always followed by a lively debate. Many students feel that the opposite 
should occur, citing the fact that, for a given perimeter P, the circle is the figure with 
the largest area. It is only when the two figures are drawn to scale, as in figure ^ that 




Figure 1. Square and circle of equal perimeter P. 



they understand the result. The point is that, for equal perimeter, the sides of the 
square lie inside the circle for most of the integration paths. 

The result can be easily generalized for any polygon with n equal sides and total 
perimeter P. figure EJ illustrates the case of n = 5. Each side has length s„ = P/n, 




x^^^ tan (n/n) 



Figure 2. Pentagon with perimeter P. (a) Pictorial representation of the vectors used 
in the calculation of A, which are defined in the text, (b) The line integral of A is 
taken along the inner (dotted) polygonal curve C„. 



placed at a distance dn = Sn/'i cot {dn/'i) from the center, where On = 27r/n. The total 
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magnetic field is simply equal to n times the field produced by a straight wire of length 
Sn carrying a current i, at a point placed at a distance (i„ from the wire, along its 
perpendicular bisector: 

^center ^ ^ j^O^ ^ J^^^2 ^-^ ^^^^y 

Substituting for = 3, 4, . . . in equation ((T)), we conclude that, for equal perimeter, 
the field at the center of a current- carrying triangle is the largest; and the fields at the 
center of other current-carrying polygons with equal perimeter decrease as the number 
of sides increases, approaching the asymptotic value of B^'^"*''^^ = -^pAn"^ obtained for 
the circle. This calculation can be assigned as a homework exercise. 

Although the area does not play a role in this example, our students usually point 
out that it should play a role in determining the auto-flux through the wire loops. For a 
given perimeter P, the areas enclosed by the polygon wires are An = cot (7r/n)/(4n), 
approaching the area of the circle, Ac = P^/(47r), as the number of sides increases. The 
naive multiplication 

p-ter^„ = /^nsin(7r/n), (2) 

grows with n. Normalizing this type of "flux" by as we shall henceforth do|, we 
find 2.6, 2.8, 3.1, and vr for n = 3, n = 4, n = 8, and the circle, respectively. This seems 
to indicate that the smaller field at the center of the circle is more than compensated 
by its larger area. Some students interpret this as a vindication of their initial intuition. 

Unfortunately, things are considerably more complicated than this simple argument 
suggests, making it interesting to revisit this problem in an advanced course on 
electromagnetism. Firstly, the magnetic field varies from point to point in space. The 
calculations of these magnetic fields may be found elsewhere for the polygon [2], for the 
circular loop jS] , and for planar wires |lj . Secondly, these fields diverge at the position of 
the wires, meaning that some regularization must be used. Thirdly, obtaining the flux 
directly from the magnetic fields requires a two dimensional integration, which becomes 
particularly difficult in the case of polygons. 

In this article, we start by calculating the vector potential A produced by a circular 
or polygonal loop of perimeter P and carrying a current i, at any point in the plane of 
the figure, inside the figure. Naturally, A and B = V x A diverge as one approaches 
the wire loop. So, we will consider the fiux of B through a surface S with edges on 
a curve C similar to (and concentric with) the current loop, but scaled down by some 
amount [c.f. figure |2fb)). Obtaining the fiux directly from B will require a further 
two-dimensional integration (besides the one needed to obtain B), which, moreover, is 
rather cumbersome in the case of polygonal surfaces. Fortunately, we may use Stokes 
theorem 

B-da= [ A-dl (3) 

s Jc 

I All our figures will be drawn for the flux $ in units of /io«P/(47r), i.e., whenever we mention $ on 
the vertical axis, we are really plotting Att^/ {fioiP). 
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to turn the two-dimensional integration involving B into the one- dimensional integration 
involving A. Many textbooks only mention the vector potential briefly; this problem 
provides a striking example of how useful the vector potential may be in practical 
applications. 

The results we obtain also provide the mutual inductance of two nested, coplanar, 
and concentric (polygonal or circular) wires of equal shape but different scales. This 
can be used for theoretical discussions and experimental studies of Faraday's law. 



2. Calculating the vector potential 

We wish to calculate An{x,y) at a point P with coordinates rp = xe^ -\- y^y, as 
illustrated in figure |2f a). We start by parametrizing the positions of the points on the 
right-hand side of the polygon as Vqi = dn^x +tey, with t G {—Sn/'2,Sn/'2). Using 
Vi = rp — Vol, we find 

4:71 ^ /■^"/2 1 droi , f"/^ dt 

— :Ani = / —dt = / , ey 

/Xo^ J-sn/2 ri dt J-s„/2 /(a; _ + {y - tf 



1 -y + + \j[x- dnf + [y- gn/2]" ^ 
-y - s„/2 + J[x- dnf +[y + Sn/2f 



The position of the points along the fc-th side (moving anti-clockwise) is simply 
given by a rotation of roi by an angle (3nk = {k — 1)6'„ = 27r(A; — l)/n. So, 
rofc = Xnk{t) ex + Ynkit) ey, where 

Xnk{t) = dn cos (3nk " t siu /3„fc, 

Ynkit) = dn sin Pnk + t COS Pnk- (5) 

As a result 

/^o^ ^J[x-X^,it)]'+[y-Y,,,it)f 

[ Sn/2-ank{x,y) + J[x-Xnk{Sn/2)f+[y-Ynk{Sn/2)]^ 

= In \ ; ^ ^ } e„fc (6) 

[ -Sn/2 - ank{x,y) + \/[x- Xnk{-Sn/2)] + [y - F„fc(-S„/2)] 

where 



e„fc = - sin Pnk Gx + cos P^k Gy (7) 



and 



±Sn/2-ank{x, y) = [x - X„fc(±s„/2)] smPnk-[y - F„fc(±s„/2)] cos/3nfc.(8) 
The final magnetic vector potential is given by 

n 

An{x, y) = J2 Ank{x, y). (9) 

k=l 
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Alternatively, we might obtain equation (jUI) from equation (0)) through the vector field 
rotations discussed by Grivich and Jackson |2]. We could now recover their Equation 
(9) with 2; = by taking B = V x A and suitable variable redefinitions§. 
As for the circular loop, we use polar coordinates. By symmetry, 

A,(p, 6) = A,{p, 6) ee = A,{p, 0) e^, (10) 

and we take rp = pe^. Parametrizing the positions of the points along the current- 
carrying circular wire of radius R a.s Vq = i?cos(/9e^ + Rsimpey, with G (0,27r), 
r = rp — Vq, and we find 



Pot Jo r cup Jo 



2^ 1 drn 



dip 



r'^'^ —R sin ipe^ + R cos ip By 
Vp^ + i?2_2pi?cosvr 



dip 



Pip + R) 



where 



K{k) 



{p' + R^)K 



dt 



'2VpR' 



{p + RfE 



'2^/pr: 



dt. 



(12) 



/o VI - kH^vT^ 

We have checked that the function An{p, 0) in equation tends to Ac{p, 0) in 
equation (jllll . as n approaches infinity. Also, by taking B = V x A and suitable 
variable redefinitions, we recover the corresponding magnetic field |3]. 



3. Calculating the flux 

We recall two points mentioned in the introduction. Because the fields diverge at the 
position of the wires, we will take the fiux in a curve similar to the original wire but 
scaled down by some amount, as in figure Efb). We may think of this as a cutoff 
introduced by the finite width of the wire, or as the situation faced in calculating the 
flux through a second loop, similar to (but smaller than) the current-carrying one. Also, 
because the direct calculation of the flux of B involves a two-dimensional integration, 
we will use equation (jS} and calculate instead the line integral of A. 

The simplicity gained in utilizing A is particularly striking in the case of the circular 
current loop, since equation (fTIH) means that A is independent of 6. Therefore, choosing 
an integration circle Cp, of radius p G (0, R), we find 



An 



An 



PoiP Jcp 



A-dl 



An 



PoiP 



A{p,0) 2np 



An 



P + R 



(p2 + R") K 



2y/pR 



{p + RfE 



2y/pR 



(13) 



P + R) \P + R, 

§ There is a subtlety concerning the fact that, since we have determined A{x,y, z) only for the plane 
z = 0, we cannot perform the derivations with respect to z. However, these do not enter the calculation 
of Bz{x, y, 0) which, by symmetry, is the only non- vanishing component of B(x, y, z) when z — Q. 



Self-magnetic flux of current- carrying polygons 



6 



where, in going to the second hne, we have made p and R dimensionless by scahng 
them by the perimeter P\\. It is instructive to compare the trivial reasoning on the first 
hne of equation with what would be needed to calculate the flux directly from the 
p- dependent Be- 

Next we consider the magnetic field produced by a polygon with perimeter P, n 
equal sides, and carrying the current i. The distance from the center to each of the sides 
is given by Consider also a second ra-sided polygon C„ whose sides lie a distance 
Xcnt £ (0, dn) from the same center. The flux through this polygon is given by 

fXcut tan (7r/n) 

y) dy.{U) 



A„.-dl 



n 



Cn 



n 



Xcut tan (-Tr/n) 



A„-dZ 

/first side 

Looking back at equation © one notices the need for integrals involving the logarithm 
of rather complicated functions. Things can be greatly simplified, however. We start 
by rescaling all distances by the perimeter P, thus rendering the variables x, y, s„, and 
dn appearing in equation © dimensionless^. Next we introduce new parameters u and 
new variables v through 

V = y~ F„fc(±s„/2), 
for use in equations (0) and Thus, for equation (jUj) we need 
Inkiu,v) = / In 



u sm 



(3nk - V COS (3nk + Vv^Tv^ 



dv - 



(15) 
(16) 



We find^ 



Ink[u,v] 



In (^—v + + + Vu"^ + 
In (v + + ? 



if Pnk 
if Pnk 







7r 



(17) 



otherwise . 



—v + u CSC (3nk In (i> + Vu'^ + 

+ {v + u cot (3nk) In (u sin f3nk — v cos j3nk + + f ^ 

Combining this with equations (jHI)^©, and substituting into equation (jHj), one obtains 



^n = n COS (3nk {ink- I, 



nk j ' 



k=l 



(18) 



where 



Ink = I [^cnt - X„fe(±S„/2), Xcut taU {jl / u) - F„fe(±S„/2)] 

-/[Xcut -X„fc(±s„/2), -Xcut tan (vr/n) - F„fe(±s„/2)] . (19) 

II We have made the variable substitutions p' — p/P and R' = R/P ^ l/(27r), and then dropped the 
primes. 

^ We have made the variable substitutions x' = x/P, y' = y/P, sj^ = Sn/P = l/n, and d'j = = 
cot (0„/2)/(2n), and then dropped the primes. 

We are very grateful to Ana C. Barroso for help with this integral. 
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We have checked equations (fT^ and (fTH|) in two important hmits. First, expanding 
around Xcnt = 0, we find that the fiuxes tend to the product of the magnetic field at the 
center with the area of a small central region whose distance to the sides is Xcut- Indeed, 

- Br'''' ir xL + 0{xl,) ^ An^ xL , (20) 
^ Br'""' n tan {t^/u) + 0{xlJ An^ tan^ (vr/n) sin (vr/n) x^ • (21) 

Here and henceforth (including in all figures), we normalize the fiuxes by fioiP / (Att) , the 
magnetic fields by fioi/ (47rP), and we continue to scale all distances by P. Naturally, we 
can recover equation ()20|) from equation ()21|) in the limit of n going to infinity. Second, 
$„ tends to $c as n goes to infinity, for all values of Xcut- This can be seen in figure 01 
which displays for n = 3, A, 8, and $c as a function of Xcut- Each fiux diverges 



18 




I — ' — , ^-^^=^-r-~r . — I — . — I — . — I — . — I — . — 1_ 
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X , 
cut 

Figure 3. Auto-fluxes as a function of a;cut, for current-carrying polygons with n = 3, 
4, 8, and for tlie circular loop. 



at 3; cut = dn, while $c diverges at Xcut = -R? providing a nontrivial crosscheck on our 
expressions. Notice that, for each value of Xcut < ds, the curve for $c lies below all other 
fiuxes. Although the fields B vary as one moves away from the center, a very rough 
way of understanding this result is the following: the field at the center Br'^'^ decreases 
as n increases — c.f. equation ((H); on the other hand, for fixed Xcut, the areas through 
which the fiux is being considered are given by nta.n{n/n) x^^^, for and by n x^^^, 
for $c, which also decrease as n increases. Therefore, in this case the "area factors" do 
not compensate for the smaller fields, as seen in equations ()20|) and (PT|) . 

Since the fiuxes diverge for Xcut = c^n, we may choose to consider another situation. 
We take all wires to be of a fixed width 6 (in units of P), and we regularize the fiuxes 
by integrating only up to p = R — 6, for the circle, and Xcut = dn — S, for the polygons. 
The results are displayed in figure^ as a function of 6, for n = 3, 4, 8, and for the circle. 
We notice the following features: i) for any finite value of S, the auto-fiux increases as n 
increases — this indicates that, here, the "area factor" is making up for the smaller value 
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of the magnetic field at the center; ii) again, the curves of tend to $c as n increases; 
iii) the flux diverges as the width of the wires tends to zero, as expected. 

Comparing figure El and figured we notice that $ decreases with n in the first case, 
while it increases with n in the second. So, in contrast to the previous case, here the 
"area factor" compensates for the smaller fields. We can get a rough understanding 
for this in the following way: for fixed 6, the areas through which the flux is being 
considered are given by 

ntan {n In) (dn — 5)^ = — cot {n In) — 5 + ntan {n In) 5^, (22) 

An 

for $„, and by 

vr(i? - 5)2 = ^ - 5 + TT (23) 
An 

for $c, in units of P^. As 5 vanishes, the areas in equations and are 
dominated by their first terms, which do increase enough as to offset the order of the 
field magnitudes. Of course, this is a very crude argument, since, because the fields 
vary in different ways as one moves away from the center, using B^^^^"^ in the reasoning 
is a considerable source of error. Nevertheless, this rough argument is consistent with 
figure m 

One can show that, although the curve of $c lies above those of for 5 7^ 0, the 
ratios $„/ $c tend to one as 5 approaches zero. This might be difficult to guess initially, 
since it seems to contradict the "area factor" , but it has an interesting interpretation in 
terms of the line integral of A. For points very close to the wires, the field approaches 
that of an infinite wire and A diverges logarithmically. Consequently, we may interpret 
the result of the line integral as the product of a logarithmic divergence with the 
perimeter P over which the integral is taken. Since these features are common to 
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all the current- carrying loops, all ratios approach unity. Of course, the same would not 
be the case if we had taken current-carrying loops of different perimeter (recall that we 
have normalized all fluxes by ^qiP / [At^)) . 

We can choose other regularizations besides the ones discussed above (constant 
Xcut and constant 5). For instance, we may ask that the surfaces through which the flux 
is being considered have the same area. In this case, as in the case of fixed Xcut, 
decreases as n increases. In contrast, if we ask that the surfaces through which the fiux 
is being considered have the same perimeter, then increases as n increases, as in the 
case of fixed 5. One can get a rough understanding for these features along the fines of 
the analysis made above. 

Finally, we recall that the line integrals of A have been performed over curves C„ 
and Cp identical to the current-carrying wires, but smaller. This is what one needs for 
the calculation of the mutual inductance between two (polygonal or circular) current- 
carrying wires of equal shape and different scales that lie on the same plane and are 
concentric. Our results apply directly to that case. 

4. Conclusions 

Motivated by a simple exercise in elementary electromagnetism, wc have studied the 
interplay between the magnetic fields and the areas of current- carrying polygonal 
and circular wires of equal perimeter. Wc have calculated the vector potential 
A for these situations, because its line integral provides a much simpler way of 
computing the magnetic fiuxes; this example illustrates the usefulness of A in practical 
calculations. Since the corresponding auto-fiuxes diverge, we have discussed a number 
of regularizations, comparing the fiuxes in each case, and seeking intuitive arguments 
for the results. As a bonus, our results can be applied directly to the calculation of 
mutual inductances in a variety of situations. 
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